A Hopf bifurcation with a strong incommensurate frequency response in quantum 

wells 
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The behavior of deha-doped wide quantum well het- 
erostructures in the presence of intense far-infrared radiation 
is studied using semiconductor Bloch equations. A quan- 
tum well is designed where one can either obtain a strong 
subharmonic (period-doubling) or a strong incommensurate 
(Hopf) frequency response by varying the sheet density and 
field strength. These strong responses are easily attainable 
with current quantum well technology and the field ampli- 
tudes and frequencies of the drive are well within the range 
of the free electron laser. 

The effective description of some quantum mechanical 
systems, far from a classical correspondence, can be non- 
linear. Strongly driven THz intersubband transitions in 
AlGaAs/GaAs quantum wells are a good example of such 
a system and have been shown to exhibit novel nonlinear 
phenomena. Many-body effects which cause the nonlin- 
earities become much more prominent in wide quantum 
wells (w 300A) where intersubband spacing is of the or- 
der of lOmeV which is approximately of the order of the 
electron-electron Coulomb interactions. The main signa- 
ture of the nonlinearities is the dynamic shielding by the 
electron gas of the incoming radiation. The screening 
blue shifts the absorption peak frequency, by an amount 
proportional to the intersubband population difference 
and sheet density, to the dressed frequency at which col- 
lective oscillations of the entire electron gas occur. This 
depolarization shift also causes the generation of second- 
harmonics of the drive frequency when it is at half of this 
dressed frequency jl]. The depolarization shift and su- 
perharmonic generation have been accurately simulated 
and measured, see In addition to this dressing 

of the intersubband frequency there is a dynamic Stark 
effect which shifts the absorption peak by an amount pro- 
portional to the intensity of the incoming radiation. Ex- 
periments by Heyman et al. ^ and Craig et al. Q showed 
the intersubband absorption peak distorting and shifting 
to the red towards the bare intersubband frequency as 
it saturates with increasing THz intensity. These results 
were in good agreement with a two-subband density ma- 
trix model proposed by Zaluzny . 

A superharmonic generation is a general feature of clas- 
sical nonlinear systems with periodic driving, whereas a 
subharmonic generation or an incommensurate frequency 
response are not as ubiquitous and usually require the 
fine tuning of parameters. The same could be the case in 



(effectively nonlinear) quantum mechanical systems and 
this necessitates a theoretical search for the appropri- 
ate parameter ranges. The generic bifurcation from a 
stationary solution or a fixed point, corresponding to a 
periodic orbit, in a classical system is a Hopf bifurca- 
tion [Q and it typically leads to a response comparable 
in magnitude to the drive. This motivates the search for 
a strong response (a Hopf bifurcation from a fixed point) 
in our effectively nonlinear quantum mechanical system 
that is a second frequency response of the collective elec- 
tron excitations incommensurate with the frequency of 
the laser drive. Theoretically this response can be ex- 
plained by a parametric oscillation of the electron gas 
where the system is able to respond at a new frequency 
that is a dressing (as explained above) of the frequency 

Another objective of this letter is to verify the pre- 
dictions of the model of Galdrikian and Birnir j2| for 
period doubling bifurcation for a more realistic situa- 
tion in the quantum well. The present model is more 
complete and makes less approximations on the electron- 
electron Coulomb interaction and also includes the ex- 
change terms in the Hamiltonian. On the other hand an- 
alytical calculations as the averaging technique used in |^ 
are considerably more difficult to perform now. There- 
fore we are limited to focus on numerical results for the 
nonlinear dynamics. 

S her win et al. [|1O|,0 have studied the optical prop- 
erties of wide quantum wells of Alo ^Gao/^As (barrier) 
and GaAs (well) driven by the far infrared radiation 
of the free electron laser (FEL) at UCSB. The epitax- 
ially grown structures are filled with a sheet density of 
sa 1.5 X lO^^e/cTO^ provided by silicon donor layers lo- 
cated outside of the well. The donors are usually at hun- 
dreds of A's away to avoid scattering of the electrons by 
the dopants, thus increasing the coherent radiation out- 
put of the electrons collective oscillations. 

Many-body effects on intersubband transitions and op- 
tical properties of QWs have been studied mainly with 
Hartree-Fock and rotating wave approximations for the 
response of two-(sub)band QWs Based on sim- 

ilar methods Nikonov, Imamoglu et al. demon- 
strated that the influence of exchange and depolarization 
terms on the absorption or gain spectrum disappears in 
the limit of very narrow quantum wells. Numerical inte- 
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gration of the multisubband semiconductor Bloch equa- 
tions were performed by Tsang et al although they 
neglected depolarization shift terms, which is a serious 
handicap for their predictions for wide quantum wells. 

In this letter we study the semiconductor Bloch equa- 
tions without RWA from a nonlinear dynamical systems 
perspective and find, in qualitative agreement with pre- 
vious work by Birnir and Galdrikian that for the 
two-subband driven near resonance period doubling bi- 
furcations or optical bistability may occur. For a three- 
subband QW designed for this project and driven at res- 
onance uj = E20 several types of bifurcations may occur 
depending on the shape of the well and the sheet density. 
We observe a Hopf bifurcation producing a strong signal 
at a lower frequency than the frequency of the drive. This 
signal corresponds to a dressing of Eiq and for a differ- 
ent (lower) value of the field strength a period-doubling 
bifurcation is observed that also produces a strong signal 
but now at half the frequency of the drive. Both signals 
are produced at moderate values of the field strength. 
The latter result is a big improvement over the result 
m |l that was much weaker than the fundamental and 
occured at relatively large field strength. Our approach 
can in principle be applied to any multisubband systems, 
but in practice when to > E20 one may run in problems 
such as coupling the system to the continuum or the ex- 
citation of LO phonons (if uj > 36meV). However, these 
problems can be avoided by introducing more barriers in 
the well which create several closely spaced subbands as 
in Figure 1. 

The optical properties of the confined electrons in the 
QW can be studied with the many-body Hamiltonian 
written in terms of field operators. The electron contri- 
bution is given by 



H,i = J d-'xip\x,t){~—h^\/^ + W{z))i,{x,t) + 

— / d^x / d'^x'ipUx,t)'tl>Ux' ,t)- j-ij^ix ,t)'ip{x,t) 

2 J J \x-x'\ 



e£{t) / d xTp {x,t)zip{x,t), 
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where W{z) is the bare QW potential. It is also assumed 
that the electric field generated from the radiation off 
the collective oscillations inside the QW are negligible 
compared with the laser light £. The field operators ex- 
pressed as a linear combination of annihilation operators 
ak,n are 

ipix^t) = 2^ — =C„(z)afc,„(t) 

k,n 

where Cn{z) is the self-consistent n-th Hartree eigenstate, 
k is the in-plane wave vector and A is the QW area. 
The Hamiltonian can be rewritten in terms of (time- 
dependent) annihilation and creation operators as 
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— ^ lJ,n,m{k)£(t)a\^^ak,m., (2) 

k,n,m 

where En{k) is the self-consistent Hartree energy, /i„_m 
is the dipole moment between the n-th and m-th states, 
/„(fc) is the Fermi occupation number for the nth- 
subband and the Coulomb interaction terms are 

fdzJ rf,/^l(^l)Cn.(^2)g„3(^2)g„,(^l)^-,|.,-.,| 
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where eg is the static dielectric constant in GaAs, which 
is approximately 13. In addition to the usual exchange 
terms that appear in the electron gas Hamiltonian we ob- 
tain depolarization shift terms in the Hamiltonian which 
are due to the electron shielding and the quasi two- 
dimensional character of the problem (it is not negligible 
in wide quantum wells with sheet densities > lO^^e/cm^). 
The interaction with the laser light is introduced through 
the electric dipole approximation which implies that any 
momentum transfer due to the THz radiation is ne- 
glected. 

Hamiltonian in hand, we proceed with the derivation 
of the equations of motion for the density matrix in the 
time-dependent Hartree-Fock approximation, which are 
also known as the semiconductor Bloch equations. They 
provide the collective optical response of the electrons to 
the laser light. In order to be more realistic we introduce 
dissipation in our equations of motion. Energy (inelastic) 
and momentum (elastic) scatterings are accounted for al- 
beit phenomenologically (by the relaxation-time approx- 
imation). 

The dynamics for the density matrix elements 
CTn,n' (k) = < al; n^k.n' > are obtained from the Heisen- 
berg equations of motion for the operator al „ap,„' fol- 
lowed by a quantum statistical average and the Hartree- 
Fock approximation. We finally obtain, see for details, 

- i^^n,n'{k) = E„^n'{k)an,n'{k) + 1 — ^ (o-„_„/ ( fc) - /n(fc)5„„/) 

J- nn' 

^ ^ Vriiri2n3ri{^)fTniri3 ~^ ^ | )cJ'ri2 n' (^) 
(g)o-nn3(fc)cr„2„4(|fc - q\) 
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(k) ~ fl^r 



(3) 



where En^n'ik) = En{k) — En'{k) and T„„/ is the exper- 
imentally measured relaxation time, see Heymann et al. 

Another approximation is to consider the effective 
masses of all subbands equal, which is fairly accurate for 
AlGaAs QW's. We also neglect the exchange terms since 
their inclusion does not modify qualitatively our results 
and the analysis becomes simpler. It can be shown ana- 
lytically that this is a good approximation for high sheet 
densities, since then the kinetic energy terms dominate 
p8[ . With these approximations and a sum over k per- 
formed, the dynamics for the averaged density matrix 
becomes 

d h 

— ifl' — an^n' ~ {En — -En')"^".™' + ^T^; (o"n,n' " l^nn') 

Ob 1 nn' 

— ''^^£it)(^flmnU,nn> — fJ-n'mO-nm) (4) 

m 

where (t„^„/ = -^J2k'^",ri'{k), Ns is the sheet density 
and crJJ^, is the equilibrium density matrix. The equa- 
tions were put in dimensionless form by dividing all terms 
by £^20- 

The QW studied was designed to have length 310 A, 
depth of 300 meV, two barrier steps of 26A in width 
and 150meV in height. Figure ^ shows the effective well- 
shape with a sheet density Ng — 3.0 x lO^^e/cm^. The 
equations (Q) for three subbands only were integrated 
using the fourth-order Runge-Kutta method with 2048 
steps per cycle of drive and the results are plotted in 
Figures 2-4. 

Figure |^ shows the bifurcation diagram of the dipole 
moment of the collective oscillations. The scaled dipole 
moment < iJ,(t) > /fiiQ is plotted as a function of the 
unitless field amplitude A = eSzio/ £20- A PDB in the 
Poincare map at A = 0.17 followed by a Hopf bifurcation 
(as the field decreases) at A = 0.378 are observed. We 
found an optimum sheet density in which the bifurcations 
occur with greatest strength and at the same time require 
small field amplitude, such information can be used as a 
guide for the experimentalist. For < 0.5 x lO^^e/cm^ 
and for Ng > 6.0 x lO^^e/cm^ the bifurcations become 
too small for practical observations. The phase portrait 
of the Poincare map in Figure ^ shows the scaled dipole 
moment plotted against its time derivative (on the ver- 
tical axis) for various values of field strength before and 
after the Hopf bifurcation. This plot is the strongest evi- 
dence that Hopf bifurcations occur for our system. Near 
A = 0.378 a supercritical Hopf bifurcation is observed as 
we decrease the field intensity, the fixed point, a periodic 
orbit with the frequency of the drive, becomes unstable 



generating a small closed orbit. Several of those orbits 
are plotted each one corresponding to a fixed value of the 
field strength, decreasing the field amplitude gives a new 
and larger closed orbit. 

The period-doubling can be understood as a two- 
photon process, see its increased strength in this let- 
ter is due to the fact that it is being assisted by the pres- 
ence of almost equally spaced (in frequency) frequencies 
-Eio and £^20- The dynamics of the electron gas corre- 
spond to a periodic orbit of twice the period of the re- 
sponse before the bifurcation. 

The large and broad peak at roughly w/2 is the main 
feature of the bifurcation as can be seen in Figure ^ which 
shows the power spectrum before and after the Hopf bi- 
furcation. The peak in Figure ^ is broad due to the differ- 
ence frequencies, dressed ±(i?2i — Eio), that are clearly 
visible on each side of the incommensurate peak. The 
half-frequency peaks (not shown) for the PDB are sharp 
in distinction By varying the sheet density and field 
strength the incommensurate peak can be tuned but this 
will be explored in Batista, Tamborenea and Birnir 
In Figure 4 we also observe a DC response that concides 
with the frequency axis. The dynamics of the electron gas 
after (decreasing A) the Hopf bifurcation corresponds to 
a quasi periodic orbit on a torus defined by the incom- 
mensurate frequency and the fundamental. At both the 
PDB and Hopf bifurcations the populations of the up- 
per subbands are suddenly increased at the expense of 
the zeroth subband. The frequency of the drive is of the 
order of 20 meV, which is in the range of the FEL (how- 
ever, E20 can be lowered by decreasing the tunneling if 
necessary). 

The realization that QWs can undergo a Hopf bifurca- 
tion (from a periodic orbit) only observed before in clas- 
sical nonlinear systems can lead to many applications in- 
cluding frequency down-converters. In our QW the lower 
subbands are almost equally spaced, but that could be 
easily changed by increasing its asymmetry and eventu- 
ally eliminating the possibility of PDBs. The asymmetry 
though cannot be too big otherwise the tunnelling is re- 
duced and we cannot guarantee that our three-subband 
system driven at resonance with E20 is really isolated 
from upper subbands and the continuum. 

The authors would like to thank Mark Sherwin, S. 
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and discussions. The idea of nonlinear bifurcations in 
quantum wells were inspired by J. D. Crawford's elegant 
exposition on bifurcation theory . The work was sup- 
ported by NSF grants DMS-9704874, DMS-0072191 and 
a DARPA-Navy grant N00014-99-10935. 
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FIG. 1. The stationary self-consistent potential for a sheet 
density of Ns — 3.0 x lO^^e/cm^ , and its first seven eigen- 
states with asymptotes set to their energies. The double bar- 
rier creates, through tunneling, three closely spaced subbands 
well isolated from upper subbands. The slight asymmetry 
enhances the nonlinear effects, but too much asymmetry re- 
duces the tunneling and with it the decoupling with upper 
subbands. 
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FIG. 2. The Poincare map of the averaged dipole 
moment in the doped QW with a sheet density of 
Ns = 3.0 X 10"e/cm^. With increasing field strength a PDB 
occurs at A = 0.17. At A = 0.378 a supercritical Hopf bi- 
furcation occurs, as the field is decreased below this value. 
The two branches of the PDB undergo a Hopf bifurcation at 
A = 0.19. 
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FIG. 3. Phase portrait view of the Hopf bifurcation in the 
quantum well. Each one of the closed orbits corresponds to 
a different value of the field amplitude, with A — 0.34 — 0.38 
and a sheet density of Nb = 3.0 x lO^^e/cm^. 
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FIG. 4. Power spectra of the dipole moment for the doped 
asymmetric quantum well near the onset of the supercritical 
Hopf bifurcation. As we decrease the field past A — 0.38 the 
low frequncy signal grows continuously. The zero frequency 
peak coincides with the y-axis and extends to zero. 
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